Peculiar Chemical Abundances in the Starburst Galaxy M82 
and Hypernova Nucleosynthesis 
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ABSTRACT 

X-ray observations have shown that the chemical abundance in the starburst galaxy M82 is 
quite rich in Si and S compared with oxygen. Such an abundance pattern cannot be explained 
with any combination of conventional Type I and II supernova yields. Also the energy to heavy 
element mass ratio of the observed hot plasma is much higher than the value resulted from normal 
supernovae. We calculate explosive nucleosynthesis in core-collapse hypernovae and show that 
the abundance pattern and the large ratio between the energy and the heavy element mass can be 
explained with the hypernova nucleosynthesis. Such hypernova explosions are expected to occur 
for stars more massive than > 20 — 25M0, and likely dominating the starburst, because the age 
after the starburst in M82 is estimated to be as short as ~ 10^ — 10'' yr. We also investigate 
pair-instability supernovae (~ 150 — SOOMq) and conclude that the energy to heavy element mass 
ratio in these supernovae is too small to explain the observation. 



Subject headings: galaxies: individual : M82 
starburst — nucleosynthesis 

1. Introduction 

M82 is the most active nearby starburst galaxy. 
Recently several exciting discoveries have been 
made for M82. X-ray observations have revealed 
the presence of intermediate mass black holes 
with masses 10^ - 10^ Mq in M82 (Matsumoto 
et al. 2001; Kaaret et al. 2001), whose loca- 
tions were found to coincide with the star clus- 
ters by SUBARU observations (Harashima et al. 
2002). Also the very energetic expanding molecu- 
lar super-bubble has been discovered (Matsushita 
et al. 2000). These findings have created lots of 
interest in the formation of black holes of vari- 
ous masses and the evolution of star burst galaxies 
(e.g., Ebisuzaki et al. 2001). 

The critically important information for under- 
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standing the evolution of star burst galaxies is the 
chemical abundances. The abundance informa- 
tion has also been provided by X-ray observations. 
Tsuru et al. (1997) observed M82 with ASCA in 
the 0.5 — 10 keV X-ray band and found that its 
spectrum can be fit with a three component model: 
a point-like hard component and extended soft- 
medium components. From the observed emission 
lines they also obtained abundances of O, Ne, Mg, 
Si, S and Fe, and found that the abundance pat- 
tern is peculiar: Si and S are much abundant than 
O and Fe compared with the solar ratio. Also the 
0/Fe ratio is almost solar. Tsuru et al. (1997) 
concluded that this abundance pattern cannot be 
reproduced with any combination of the previous 
Type labc and Type II supernova (SN labc and 
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SN II) yields (Nomoto, Thielemann & Yokoi 1984; 
Woosley & Weaver 1995; Thielemann, Nomoto & 
Hashimoto 1996; Nomoto et al. 1997a,b). They 
also discussed that the energy to heavy clement 
mass ratio of the observed hot plasma is too large 
to be of the supernova origin. 

We should note that Tsuru et al. (1997) as- 
sumed that all stars above 8Mq have already ex- 
ploded. However, the age of the starburst is esti- 
mated to be ^ lO^yr from the radio observation 
(Matsushita et al. 2000) and ~ 10'' yr from the size 
of the X-ray halo observed with ROSAT (Strick- 
land, Ponman & Stevens 1997). Only the very 
massive stars have exploded in such a short time. 

Moreover, the previous SN II and SN Ibc yields 
adopted in Tsuru et al. (1997) were obtained only 
for the explosion energy of 10^^ erg. Recently 
some massive supernovae have been found to ex- 
plode much more energetically as "hypernovae" 
than normal SNe II (e.g., Iwamoto et al. 1998; 
Nomoto et al. 2001). Therefore, we need to re- 
consider the previous conclusions by Tsuru et al. 
(1997). 

In this paper we re-examine whether the abun- 
dance and energetics of M82 are consistent with 
the supernova models taking account of the star- 
burst age and hypernova explosions. Because of 
the short age, the abundance is likely to be de- 
termined mainly by relatively massive supernovae. 
The contributions of high energy explosions would 
affect the energy to the heavy element mass ratio 
and the abundance pattern in the galactic winds. 
Recently we have found that in nucleosynthesis of 
hypernovae, the Si and S abundances are much 
enhanced relative to O; and also large Fe/0 ratio 
([Fe/0] ~ 0) can be realized (Umeda, Nomoto & 
Nakamura 2000; Nakamura et al. 2001). These 
patterns are consistent with those observed in 
M82, which has motivated us to calculate detailed 
nucleosynthesis in massive energetic core-collapse 
SNe for comparison with the M82 data. We also 
investigate pair instability supernovae (PISNe) of 
150 — 3OOM0 stars, because PISNe also yields rel- 
atively abundant Si and S compared with O, and 
large amount of Fe. 

2. Abundance Pattern 

The heavy element abundances in M82 ob- 
served with ASCA (Tsuru et al. 1997) are shown 



in Figure 1 by open circles with error bars. For Ar, 
only the upper limit is given. Here, the heavy ele- 
ment abundances are defined by the mass fractions 
of elements with respect to light elements (mostly 
hydrogen and helium). The data are normalized 
to the solar ratios (Anders & Grevesse 1989). We 
note two features of the observational data. One 
is the overabundances of Si and S relative to O 
with respect to the solar ratio, i.e., [Si,S/0] ~ 0.9. 
Here, [X/Y] = log(X/Y)- log (X/Y) q. Second 
is such low abundances of O and Fe as [0/H]~ 
[Fe/H] ~ -1.2. 

Since these abundances were derived by fit- 
ting the ASCA spectra with the two-temperature 
plasma model for the hot gas in M82, one might 
argue that these features are an artifact due to the 
too simple model assumption. Therefore, we fit- 
ted the ASCA spectrum with a multi-temperature 
plasma model in which an emission measure of 
each temperature component follows a power-law 
in temperature (i.e., the emission measure from 
the temperature T is proportional to (T/Tmax)", 
where Tmax and a are free parameters). We ob- 
tained a = 0.20 ±0.36 and T^ax = 0.93 ±0.14 keV 
as the best-fit results, and the best-fit abundances 
are quite similar to those of the two-temperature 
model fitting. We thus conclude that the above 
two features are not the artifact. 

To compare with the observed abundance pat- 
tern of M82, we calculate nucleosynthesis in the 
SN II explosions of 2OM0 stars for various ener- 
gies and metallicities. (More details on our mas- 
sive star evolution, explosion simulation and nu- 
cleosynthesis calculations are given in Umeda et al. 
2000 and Umeda & Nomoto 2002). Since [Fe/H] 
^ —1, it is better to use low metal star models to 
compare with the M82 abundances. 

2.1. Normal Type II Supernovae 

First, let us compare the abundance pattern 
of M82 with the normal SN II models for which 
the explosion energy is assumed to he E51 = 
E/lO^^evg = 1. This is the typical energy for SNe 
II as estimated from SNe 1987A, 1993J, and 19941 
(e.g., Shigeyama & Nomoto 1990; Blinnikov et al. 
2000; Nomoto et al. 1993, 1994). The ejected 
mass of Fe (mostly ^^Fe produced by the ''^Ni de- 
cay, and some ^'*Fe and ^^Fe) depends on the un- 
certain "mass-cut", but here we assume that the 
ejected ^^Ni mass is O.O7M0 as estimated for the 
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above SNc. The mass-cut depends on the progen- 
itor mass and the explosion energy. 

In Table 1, we summarize the ejected masses 
of selected elements for a various set of metal- 
licity Z and the mass M of SN II models with 
£^51 = 1. The abundances in the ejecta depend on 
the convective efficiency and the -^^C(a,7)^^0 
rate adopted for the stellar evolution models as 
shown in Table 1. Here, the parameter is pro- 
portional to the diffusion coefficient of convective 
mixing (Spruit 1992; Umeda et al. 1999). In Fig- 
ure 1, we plot the abundances of these 20Mq SNe 
II models for various metallicities. The left and 
right panels show the cases with relatively slow 
and fast convective mixing (/j,=0.05 and 0.15), 
respectively. The overall abundances are normal- 
ized to the observed Si abundance. In this section, 
we can compare only the theoretical and observed 
abundance ratios. The overall normalization is 
constrained by the energy to heavy element mass 
ratios discussed in the section 3. 

Prom Figure 1, we find that these models are 
all have much smaller Si/0 ratios than the M82 
data. The Si/0 ratios are larger for = 0.05 than 
fk =0.15 but still smaller than the observation by 
a factor of ^4. For other £^51 = 1 models in Table 
1, we find that [Si/0] < 0.57, being much smaller 
than the observation in all models. Generally, 
more massive core-collapse SNe (M < I3OM0) 
yield smaller Si/0 ratios because the O yield in- 
creases more rapidly with increasing mass than the 
Si yield. Therefore, the Si/0 ratio is likely to be 
largest for the smallest mass SN II. However, as 
shown in Table 1, the Si/0 ratios for the 13M0 
SN II models arc not large enough to be consistent 
with the M82 data. Contribution of even smaller 
mass SNe II would be negligible unless the Initial 
Mass Fimction was extremely steep. Stars with 
M > laOAf© explode as PISNe, and the nucle- 
osynthesis pattern may be different. These cases 
are discussed in Section 4. This result is quite gen- 
eral, being consistent with our previous and other 
groups' results. We also note that the Si/0 ratio 
is insensitive to the ^^C(a, 7)^^0 rate (Nomoto et 
al. 1997a). 

It is interesting to note that, the abundance ra- 
tios for these even Z elements are insensitive to the 
metallicity (Figure 1). This is because the synthe- 
sis of these elements are mostly governed by the a 
- nuclei reactions, which are almost independent of 



metallicity (Umeda et al. 2000). From this reason, 
in the rest of this paper, we use the Z=0 models 
for the progenitor and supernova models. 

2.2. Hypernovae 

Recent observations suggest that some massive 
SNe with M > 20 - 25M0 explode more than ten 
times energetically than normal SNc II and SNe 
Ibc, which may be called "hypernovae" . The first 
discovery of such SNe was SN1998bw (Galama et 
al. 1998; Iwamoto et al. 1998). Then several 
other hypernovae and their candidates have been 
discovered, namely, SNe 1997ef, 1997cy, 1999as, 
and 2002ap (e.g., Nomoto et al. 2001; Mazzah et 
al. 2002). 

Umeda et al. (2000) and Nakamura et al. 
(2001) have shown that the large Si/0 ratio may 
be the signature of the energetic supernova explo- 
sions. Figure 2 compares the post-explosion abun- 
dance distribution of the hypernova and the nor- 
mal SN II models. In the more energetic explosion, 
the outer boundary of the explosive O burning re- 
gion moves outwards. The explosive O burning to 
produce Si, S and Ar takes place when the peak 
temperature after the supernova shock passage ex- 
ceeds Tg ~ 4, where Tg is the peak temperature in 
units of lO^K. The region after the shock passage 
is radiation dominant, so the peak temperature is 
approximately related to the stellar radius r and 
the explosion energy £ as £" ~ aT'^Anr^ /3, i.e., 
Tg = El(^{r/3.16 X 10* km)-^/*. For larger E, 
therefore, the O burning region extends to larger 
r. On the other hand, the outer edge of the 0-rich 
region is fixed at the C+O/He interface, where the 
density is too low for explosive burning. This is 
why the Si/0 ratio is enhanced for more energetic 
explosions. 

Figure 3 shows the nucleosynthesis pattern of 

the 25M0 {Z=0, fk = 0.05) models for vari- 
ous explosion energies, £51=1, 10 and 30. Here 
the abundances are normalized to the observed Si 
abundance and the mass-cuts (i.e., the compact 
remnants' masses) are chosen to eject the Fe mass 
to fit the data point. Thus the ejected Fe masses 
are 0.07, 0.095 and 0.12 Mq, for £51=1, 10 and 30, 
respectively. The mass-cuts for these models are 

2.3, 3.1 and 4.0 Mq, respectively. Table 2 sum- 
marize the ejected masses of O, Ne, Mg, Si, S and 
Fe as a function of the explosion energy for the 25 
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and 30 M© {Z=0) models. 

The observational and theoretical mass ratios 
of Si/0 are compared in the top panel of Figure 4. 
The observed large abundance ratio of Si/0 can 
be reproduced for £^51 > 10 in the 25Mq models, 
and E^i > 30 in the 3OM0 models. Here, the 
results are for the fk = 0.05 models. For the fk = 
0.15 models, which is more 0-rich, higher energy 
is required to be consistent with observations. 

Similar Si/0 ratios can be obtained for more 
massive stars if the explosion energy is sufficiently 
large. Therefore, more massive core-collapse hy- 
pernovae with initial masses of ^ 50 — lOOM© 
could also be consistent with the M82 abundance 
pattern, although we do not know how much frac- 
tion of these massive stars actually explode rather 
than collapsing to black holes without explosion. 

In these results, Ne and Mg appear to be un- 
derproduced. However, the ratios (Ne, Mg) / O are 
sensitive to the uncertain reaction rate ^'^C{a, 7)"'^^0 
and convective parameter fk as seen in Figure 1. 
Thus a better agreement with observations can 
be obtained by choosing appropriate parameters, 
which does not change the present conclusion on 
the Si/0 ratio. 

3. Energy to Heavy Element Mass Ratio 

In this section, we discuss another observational 
constraint which theoretical models should satisfy. 
Suppose that the hot plasma in M82 is a mix- 
ture of the SN ejecta and the interstellar gas, and 
most of the thermal energy of the hot plasma was 
supplied by supernovae. Then we can constrain 
the ratio of the energy released by one supernova, 
E, to the mass of the heavy elements as follows 
(Tsuru et al. 1997). 

The thermal energy of the hot plasma supplied 
by one SN explosion can be written as 

2 fj/mu 

where M^j and Mam are the masses of the ejecta 
and the ambient gas heated by a single SN and, 
0.6) is the mean molecular weight. The ratio 
of the mass of each heavy element (M*) to the 
total mass (Mtotai) is expressed as 

M* M* + M* M* 

-1" ej ' ^"am ^ -'"ej /2-1 



where the quantities with asterisk denote the mass 
of each heavy element. From these equations, we 
obtain: 

^ > 4 9x 10-3 (^] (^] " = ( ^] 

(3) 

where M*^q is the ejected mass in units of Mq of 
each heavy element synthesized by the SN. For a 
given heavy element mass in the ejecta, this in- 
equality determines the minimum explosion en- 
ergy required for the supernova model to be con- 
sistent with the observations. 

Using this relation and the observed values of 
kBT(~ 1 keV) and M*/Mtotid (Tsuru et al. 1997 
and references therein), the minimum energy to 
the heavy element mass ratios for O, Si, and Fe 
are given as follows: 

(£5i/M;jQ)mi„(0) = 7.9tli, (4) 
(£5i/M*Q)mi„(Si) = 17.5tl;^, (5) 

(i?5i/M*Q)mi„(Fc) ^ 51.2tlli. (6) 

The limits for other elements can also be obtained, 
but these three elements are sufficient for con- 
straining E^i. 

Here, we have neglected the contribution of the 
stellar wind to the total energy of the hot plas- 
mas. The energy of stellar wind from solar metal- 
licity stars is estimated to be ~ 0.2 x 10^^ erg 
from the observation of OB associations (Abbott 
1982). This is much smaller than the typical ex- 
plosion energy of a normal SN II, ~ 10^^ erg. This 
energy is even smaller for a lower metal environ- 
ment as in M82, because the wind mass-loss rate 
decreases with metallicity as oc (Leitherer, 
Robert & Drissen 1992). Furthermore, we are con- 
sidering the energy from a hypernova with > 10^^ 
erg. Therefore, it is not a bad approximation to 
neglect the energy from stellar winds in the hot- 
plasma. 

In the middle and bottom panels of Figure 4, we 
compare the observed ratios (i?5i/M*jQ)niin with 
the theoretical E^i/M*-^q. For the Er,i = 1 mod- 
els, the ratios £^51 /M*-q for all elements are much 
smaller than the minimum values in Eq. (3). On 
the other hand, all other energetic models shown 
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here satisfy the constraints. Therefore, the 25Mq 
models with E51 > 10 and the SOM© models with 
-£^51 ^30 are consistent with both constraints on 
the Si/0 abundance ratio and the energy to mass 
ratio. Here, we note again that these results are 
for the fk = 0.05 models and larger energy is re- 
quired for fk = 0.15 models to be consistent with 
observations. 

4. Pair Instability Supernovae 

It has been known that the stars with M ~ 
150 — 3OO-M0 also explode energetically as pair 
instability SNe (PISNe) (e.g., Barkat, Rakavy & 
Sack 1967; Fraley 1968; Ober, El Eid & Pricke 
1983; Woosley & Weaver 1982). During oxygen 
burning the cores of these stars enter into the 
electron - positron pair creation region. Then 
the adiabatic index of these cores becomes less 
than 4/3 and the stars collapse rapidly. The cen- 
tral temperature increases during the collapse and 
rapid 0-burning takes place to produce a large 
amount of thermal energy and Si. For the stars 
with M < 3OOM0 the thermal energy produced 
by oxygen burning exceeds the gravitational po- 
tential energy, and the stars disrupt completely. 
Compared with normal SNc II, the explosion en- 
ergies are larger, typically 2 — 4 x 10^^ erg, and 
the Si abundance in the ejecta is also larger (see 
Table 2). During the PISN explosion a large 
amount of Fe (~ IOM0) can be produced (Umeda 
& Nomoto 2002; Heger & Woosley 2002), although 
the amount of Fe production is very sensitive to 
the maximum temperature attained. 

In Figure 5, we show the abundance distri- 
bution of the 2OOM0 PISN model. In Table 2, 
the ejected masses of O, Ne, Mg, Si, S and Fe 
are shown for the 170 and 2OOM0 Pop III PISN 
(Umeda & Nomoto 2002). The abundance pat- 
terns of the 170 and 2OOM0 stars are compared 
with the M82 data in Figure 6. It seems that the 
large [Si/0] ratio in M82 can be reaUzed by PISN 
nucleosynthesis . 

However, we also need to compare the energy to 
mass ratio for the PISN models. Since the amount 
of ejected O is least model dependent, we show 
the energy to the O mass ratio for the PISNe. 
Our models with M= 170 and 2OOM0 {Z=0) have 
E5i/M*.q{0)= 22.4/44.2 = 0.51 and 26.8/56.0 = 
0.48, respectively. These are much smaller than 



the observed lower limit of 7.9 (Eq. 4), because 
the ejected O mass is more than ten times larger 
than in the core-collapse SNe. 

This ratio may be larger for more extremely 
massive stars. For example, in the model of Heger 
& Woosley (2002), their I3OM0 He star model, 
corresponding to the initial stellar mass of 26OM0, 
gives E5i/M*.q{0)= 87.1/29.9 = 2.9. For these 
models, however, the O /Fe ratios are typically too 
small ([0/Fe]=-l.l for the I3OM0 He star model) 
to be compatible with observations. 

5. Conclusions cind Discussion 

X-ray observations have shown that the chemi- 
cal abundance in the starburst galaxy M82 is quite 
rich in Si and S compared with oxygen. Such 
an abundance pattern cannot be explained with 
any combination of conventional Type I and II 
supernova yields. Also the energy to heavy el- 
ement mass ratio of the observed hot plasma is 
much higher than that of normal supernovae. We 
have calculated explosive nucleosynthesis in core- 
collapse hypernovae to show that the abundance 
pattern and the large energy of the hot plasma 
of M82 can be explained with such hypernovae as 
(M, ^51) = (25, 10) and (30, 30). More mas- 
sive and more energetic core-collapse explosions 
would also satisfy the observed constraints. We 
have also investigated pair-instability supernovae 
(M ~ 150 — 3OOM0) and conclude that the energy 
to heavy element mass ratio in these supernovae 
is too small to explain the observation. 

Such "hypernova" (energetic core-collapse SN) 
explosions are expected to occur for stars more 
massive than > 20 - 25M0 (MazzaU et al. 2002). 
The upper mass limit of core-collapse SNe is still 
uncertain. Too massive stars may collapse to black 
holes without explosion. 

The question is how the abundance in M82 can 
be dominated by such hypernovae. One possible 
explanation is the age effect. As mentioned in In- 
troduction, the age after the beginning of star- 
burst is estimated to be ~ 10^ — 10^ years. On 
the other hand, the lifetime of our Pop HI stars 
are 1.40, 1.10, 0.81 and 0.70 xlO^ years for 15, 
20, 25 and 3OM0 models, respectively (Umeda et 
al. 2000). This is consistent with the assumption 
that only such massive stars as M > 20 — 25M0 
have exploded in M82. 
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There have been some suggestions that hyper- 
novae might be correlated with gamma-ray bursts 
(GRBs) (e.g., Iwamoto et al. 1998). However, how 
much fraction of hypernovae is associated with 
GRBs is still unknown. Paczynski (2001) dis- 
cussed that the SN rate (including the energetic 
SNe) is much higher than the GRB rate. On the 
other hand, we have shown that in order to explain 
the observed abundance pattern of M82, dominant 
fractions of massive stars (M >^ 25Mq) needs 
to be hypernovae. This implies that hypernovae 
are much more frequent than GRBs. If GRBs 
are associated with hypernovae, it is likely that 
they occur only for certain special cases of hyper- 
nova explosions. For example, the explosion en- 
ergy needs to exceed a certain value, the explosion 
needs to be extremely aspherical, i.e., the ejecta 
needs to be strongly beamed, and both hydrogen 
and helium envelopes need to be stripped off be- 
fore the explosion. In fact, GRB980425 associate 
with SN1998bw is the exceptionally weak GRB 
(Galama et al. 1998), while SN1998bw is among 
the most energetic hypernovae so far (Nomoto et 
al. 2001; MazzaU et al. 2002). 

What is the implication of this work to the in- 
termediate mass black holes recently discovered in 
M82? One possible scenario for the formation of 
such black holes is the merging of massive stars in 
the dense stellar clusters (Portegies et al. 1999). 
This scenario is consistent with the results of this 
paper, because such merging increases the num- 
ber of massive star explosions. However, our re- 
sults constrain the merging history. The typical 
mass range for the merged stars cannot be in the 
range for PISN (M ~ 150 - 300Mq) to be consis- 
tent with the energy to heavy element mass ratios 
observed in M82. 

This work has been supported in part by the 
grant-in-Aid for Scientific Research (12640233, 
14047206, 14540223) of the Ministry of Educa- 
tion, Science, Culture, Sports, and Technology in 
Japan. 
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Fig. 1. — Abundance patterns in the ejecta of 
20Mq SN II models for various metallicity com- 
pared with the observed abundances (relative to 
the solar values) of M82. Here the left and 
right panels show the cases with the relatively 
slow convective mixing (//c=0.05) and fast mix- 
ing (/fe=0.15), respectively. The open circles with 
error bars show the M82 data. Theoretical abun- 
dances are normalized to the observed Si abun- 
dance. 
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Fig. 2. — Abundance distribution after the super- 
nova explosion of normal (E^i—l) and energetic 
(i?5i=10) core-collapse SNe {fk = 0.05 models). 
The unstable isotope '^^Cr decays into ''^Ti. 
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Fig. 3. — Same as Figure 1 but for 25 and 
30Mq Pop III SN II with normal explosion en- 
ergy i?5i = 1 and hypernova models with E^i = 
10, 30 (25Mo) and £51^ 30, 50 (3OM0). 
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Fig. 4. — Comparison of the theoretical Si to 
O mass fraction ratios with the observed range, 
which is shown by the arrows, as a function of 
explosion energy (top panel). Comparison of the 
theoretical energy to heavy element (Fe, Si and 
O) mass ratios with the observed lower limits as 
a function of the explosion energy (middle panel 
for the 25Mq models and bottom panel for the 
3OM0 models, respectively). Here, shown are the 
/fe — 0.05 models. For the fk = 0.15 models, 
which are more 0-rich, higher energy is required 
to be consistent with observations. 
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Fig. 5. Same as Figure 2 but for 2OOM0 Pop III 
PISN model. 




Fig. 6. Same as Figure 1 but for 170 and 2OOM0 
Pop III PISN models. 



Table 1 

Ejected masses of selected elements in Mq for various metallicity Z and mass M in SN II 

MODELS WITH EXPLOSION ENERGIES E^i = 1. 



Ejected Mass (Mq) 



(M, Z) 


C 


N 


O 


Ne 


Mg 


Si 


S 


Ar 


Ca 


Fe 


[Si/O] 


(13, 0.02) 


0.11 


0.058 


0.26 


0.037 


0.025 


0.072 


0.039 


0.0074 


0.0059 


0.091 


0.57 


(13, 0) 


0.19 


0.10 


0.22 


0.026 


0.023 


0.049 


0.026 


0.0051 


0.0048 


0.070 


0.56 


(15, 0.02) 


0.16 


0.060 


0.39 


0.058 


0.035 


0.085 


0.044 


0.0081 


0.0066 


0.094 


0.47 


(15, 0) 


0.19 


0.014 


0.53 


0.072 


0.030 


0.11 


0.063 


0.013 


0.013 


0.070 


0.43 


(20, 0.01) 


0.31 


0.028 


1.76 


0.12 


0.076 


0.28 


0.15 


0.026 


0.022 


0.082 


0.33 


(20, 0.004) 


0.33 


0.015 


1.27 


0.10 


0.064 


0.23 


0.13 


0.024 


0.021 


0.081 


0.39 


(20, 0.001) 


0.29 


0.0034 


1.65 


0.14 


0.077 


0.28 


0.17 


0.032 


0.030 


0.073 


0.36 


(20, 10""*) 


0.26 


4.8e-4 


1.27 


0.18 


0.080 


0.28 


0.16 


0.028 


0.024 


0.071 


0.48 


(20, 10"^) 


0.26 


0.0018 


2.02 


0.15 


0.098 


0.33 


0.21 


0.041 


0.040 


0.070 


0.35 


(20, 0) 


0.26 


2.7e-4 


1.56 


0.12 


0.071 


0.25 


0.15 


0.030 


0.030 


0.070 


0.34 


(25, 0.02) 


1.01 


0.050 


5.79 


0.66 


0.22 


0.67 


0.30 


0.049 


0.032 


0.11 


0.20 


(25, 0) 


0.61 


3.6e-4 


2.18 


0.20 


0.12 


0.31 


0.17 


0.033 


0.031 


0.070 


0.28 


(30, 0.02) 


0.96 


0.061 


3.43 


2.31 


0.68 


0.25 


0.091 


0.017 


0.0083 


0.12 


-0.01 


(30, 0) 


0.36 


1.9e-4 


4.80 


0.27 


0.24 


0.49 


0.27 


0.052 


0.050 


0.070 


0.14 


(40, 0.02) 


0.81 


0.086 


6.41 


2.06 


0.48 


0.36 


0.16 


0.026 


0.019 


0.12 


-0.12 


(20, 0.02)» 


0.27 


0.048 


2.03 


0.65 


0.17 


0.15 


0.068 


0.011 


0.0084 


0.10 


-0.01 


(20, 0)" 


0.21 


9.9e-4 


1.77 


0.47 


0.11 


0.16 


0.089 


0.018 


0.017 


0.071 


0.09 


(50, 10-*)'' 


5.89 


0.0015 


11.6 


0.38 


0.21 


0.37 


0.17 


0.030 


0.027 


0.070 


-0.36 


(50, 10-'')<= 


1.14 


9.7e-4 


12.5 


2.41 


0.79 


0.95 


0.45 


0.077 


0.068 


0.071 


0.01 



Note. — The mass of Fe (^'^'Fe, ^®Fe, ^^Fe) depends on the mass-cut. Here, the mass-cut is chosen to eject 
0.07 Mq of ^®Ni, which decays into ^^Fe, via Co. Except for the last four items, the adopted parameters are 
/fc=0.05 and 1.4 X Caughlan & Fowler (1988: CF88 hereafter) for the ^^C(a,7)^''0 rate. 

''/fc=0.15, 1.3 X CF88 for the ^^C(a,7)^'^0 rate 

''/ji=0.1, CF88 i^C(a,7)i''0 rate 

■=^=0.5, CF88 ^^C(a,7)^''0 rate 



Table 2 

Same as Table 1 but for the Z=0, M=25 & 30 Mq core-collapse SNe with various 

EXPLOSION energies. 



Ejected Mass (Mq) 



(M, Esi) 


C 


N 


O 


Ne 


Mg 


Si 


S 


Ar 


Ca 


Fe 


[Si/O] 


(25, 1) 


0.61 


3.6e-4 


2.18 


0.20 


0.12 


0.31 


0.17 


0.033 


0.031 


0.07 


0.28 


(25, 10) 


0.40 


7.2e-4 


1.55 


0.11 


0.13 


0.43 


0.24 


0.048 


0.044 


0.095 


0.58 


(25, 30) 


0.17 


0.0023 


1.03 


0.059 


0.12 


0.54 


0.30 


0.060 


0.052 


0.12 


0.85 


(30, 1) 


0.36 


1.9e-4 


4.80 


0.27 


0.24 


0.52 


0.27 


0.052 


0.050 


0.11 


0.17 


(30, 20) 


0.24 


2.2e-4 


3.21 


0.14 


0.15 


0.75 


0.44 


0.079 


0.074 


0.16 


0.50 


(30, 30) 


0.17 


1.4e-4 


2.77 


0.10 


0.13 


0.78 


0.51 


0.093 


0.081 


0.16 


0.58 


(30, 50) 


0.11 


9.1e-5 


2.16 


0.058 


0.086 


0.78 


0.56 


0.13 


0.092 


0.18 


0.69 


(170, 22.4) 


2.30 


0.010 


44.2 


1.19 


1.94 


16.2 


8.06 


1.42 


1.32 


3.6 


0.69 


(200, 26.8) 


4.24 


5.8e-4 


56.0 


3.75 


3.08 


21.2 


13.1 


2.36 


2.32 


7.2 


0.71 



Note. — Here, the mass-cut (or Fc mass) is chosen to satisfy the observed Si/Fc ratio. The M— 170 & 200 
Mq PISN models are also shown. The SN II models with E^i = 1 and PISN models do not satisfy the 
observational constraints Eq.s (4-6), but the core-collapse hypernova models (E51 > 10) shown here satisfy 
the constraints. The adopted parameters are /fc=0.05 and 1.4 x CF88 for the i^C(a, 7)1^0 rate. 
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